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ABSTRACT 

Interferometric observations of the W33A massive star-formation region, performed 
with the Submillimeter Array (SMA) and the Very Large Array (VLA) at resolutions 
from 5" (0.1 pc) to 0.5" (0.01 pc) are presented. Our three main findings are: (1) 
parsec-scale, filamentary structures of cold molecular gas are detected. Two filaments 
at different velocities intersect in the zone where the star formation is occurring. This is 
consistent with triggering of the star-formation activity by the convergence of such fila- 
ments, as predicted by numerical simulations of star formation initiated by converging 
flows. (2) The two dusty cores (MM1 and MM2) at the intersection of the filaments are 
found to be at different evolutionary stages, and each of them is resolved into multiple 
condensations. MM1 and MM2 have markedly different temperatures, continuum spec- 
tral indices, molecular-line spectra, and masses of both stars and gas. (3) The dynamics 
of the "hot-core" MM1 indicates the presence of a rotating disk in its center (MM1- 
Main) around a faint free-free source. The stellar mass is estimated to be ~ 10 M@. A 
massive molecular outflow is observed along the rotation axis of the disk. 

Subject headings: H II regions - ISM: individual objects (W33A) — stars: formation 



1. Introduction 

Stars form by accretion of gas in dense molecular-cloud cores. However, the differences, if any, 
in the details of the formation process of massive stars (those with roughly M* > 8 Mq) compared 
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to low -mass stars are not well underst ood. Recent reviews on the topic are those bv lBeuther et al 
(j2007h . and lzinnecker fe Yorkel (|2007h . 



We are carrying out a program aimed at studying how the formation of massive stars in clusters 
proceeds in the presence of different levels of ionization, from the onset of detectable free-free 
emission to the presence of several bright ultracompact (UC) H II regions. In this paper we present 
our first results on the mass ive star-formatio n region W33A (also known as G12.91-0.26), at a 
kinematic distance of 3.8 kpc (jJaffe et al.lll982l ). W33A is part of the W33 giant H II region complex 
(j Wester houtl 1 1 958l ) . It was recognized as a region wit h very high far in fr ared luminosity ( ~ 1 x 10 5 



Lq ), but very faint radio-continuum emission by IStier et al,l ()1984l ). Ivan der Tak et al.l (|2000l ) 
modeled the large-scale (arcminute) cloud as a spherical envelope with a power-law density gradient, 
based on single-dish mm/submm observations. Those authors also presented mm interferometric 
observations at several-arcsecond resolution that resolved the central region into two dusty cores 
separated by ~ 20, 000 A U. The brightest mm cor e contains faint (~ 1 mJy at cm wavelengths) 
radio-continuum emission (IRengarajan fe Holll996l ) resolve d at 7 mm into possibly thre e sources 
separated by less than 1" (rs 40 00 AU) from each othe r (Ivan der Tak Sz Mentenl 120051 ) . These 



radio sources were int erpret e d by Ivan der Tak fc Mentenl JiooJ as th e gravitationally t rapped H 



II regions set forth by iKetol (|2003l ). However, the earlier detection by iBunn et al.l (|1995) of near- 
infrared recombination line (Bra) emission with FWHM = 155 km s _1 suggests that at least some 



of the radio free-free emission is produced by a fast ionized outflow. More recently, iDavies et al 



(|2010l ) reported spectroastrometry observations of Br7 emission toward W33A. The Br7 emission 
appears to be produced by at least two physical components: broad line wings extending to a few 
hundreds of kilometers per second from the systemic velocity appear to trace a bipola r jet on scales 



of a f ew AU, while the narrow-line emission may be attributed to a dense H II region (jDavies et al 



2010) . Being a bright mid- and far-infrared source, W33A has also been target of interferometry 
experiments at these wavelengths, which reveal density gradients and non -spherical geometry in 
the warm dust within the inner few hundred AU (|de Wit et al.l 120071 . l2010l ). 



Here we report on millimeter and centimeter interferometric observations performed with the 
Submillimeter Array (SMA) and the Very Large Array (VLA) at angular resolutions from ~ 5" to 
0.5". We find a massive star-forming cluster embedded in a parsec-scale filamentary structure of 
cold molecular gas. The dense gas is hierarchically fragmented into two main dusty cores, each of 
them resolved into more peaks at our highest angular resolution. The main cores appear to be at 
different evolutionary stages, as evidenced from their differing spectra, masses, temperatures, and 
continuum spectral indices. The warmer core harbors faint free-free emission centered on a rotating 
disk traced by warm molecular gas. The disk powers a massive molecular outflow, indicating active 
accretion. In Section 2 of this paper, we describe the observational setup. In Section 3 we list our 
results, in Section 4 we present a discussion of our findings, and in Section 5 we give our conclusions. 
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Observations 
2.1. SMA 



We observed the W33A region with the Submillimeter Arraj0 (|Ho et al. 2004 ) in the 1.3-mm 



(230 GHz) band using two different array configurations. Compact-array observations were taken 
on 2007 July 17, and covered baselines with lengths between 7 and 100 KA (detecting spatial 
structures in the range of 29.5" to 2.1"). Very Extended (VEX) configuration data were taken on 
2008 August 2, with baseline lengths from 23 to 391 KA (9.0" to 0.5"). For both observations, the 
two sidebands covered the frequency ranges of 219.3 — 221.3 and 229.3 — 231.3 GHz with a uniform 
spectral resolution of ~ 0.5 km s . 

We also report on the continuum emission from archival observations taken in the 0.9-mm (336 
GHz) band on 2006 May 22. The array was in its Extended configuration, with baseline lengths 
from 18 to 232 KA (11.4" to 0.8"). These data were used to constrain the spectral index of the 
continuum sources. 

The visibilities of each data set were separately calibrated using the SMA's data calibration 
program, MIR. We used Callisto to obtain the absolute amplitude and quasars to derive the time- 
dependent phase corrections and frequency-dependent bandpass corrections. Table 1 lists relevant 
information on the calibrators. We estimate our flux-scale uncertainty to be better than 15%. 
Further imaging and processing was done in MIRIAD and AIPS. 

The continuum was constructed in the (u, v) domain from the line-free channels. The line- 
free continuum in the 1.3-mm Compact-configuration data was bright enough to perform phase 
self-calibration. The derived gain corrections were applied to the respective line data. No self- 
calibration was done for the higher angular resolution data sets. 



2.2. VLA 

We observed the (J,K) = (1,1) and (2,2) inversion transitions of NH3 with the Very Large 
ArrajJE Observations were carried out on 2004 June 14 and 15 (project AC733). The array was in 
its D configuration, with baseline lengths in the range of 3 to 79 KX (detecting scales from 68.7"to 
2.6"). The correlator was set to the 4-IF mode. Each of the IF pairs was tuned to the (1,1) and 
(2,2) lines, respectively, covering a bandwidth of 3.1 MHz (39 km s -1 ) at a spectral resolution of 
0.6 km s _1 . 

1 The Submillimeter Array is a joint project between the Smithsonian Astrophysical Observatory and the Academia 
Sinica Institute of Astronomy and Astrophysics and is funded by the Smithsonian Institution and the Academia Sinica. 

2 The National Radio Astronomy Observatory is operated by Associated Universities, Inc., under cooperative 
agreement with the National Science Foundation. 
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The data were calibrated and imaged using standard procedures in the AIPS software. Table 
1 lists the quasars used to derive the absolute flux scale, the time-dependent gain corrections, and 
the frequency-dependent passband calibration. The absolute flux scale is accurate within a few 
percent. No self-calibration was performed. 



3. Results 
3.1. Continuum Emission 

3.1.1. Morphology 



Our observations at 1.3 mm resolve each of the two mm cores reported by Ivan der Tak et al 



(|2000j) into multiple continuum sources. The concatenation of Compact- and VEX-configuration 
data permits us to simultaneously resolve the structures at ~ 0.5" resolution and to be sensitive to 
relatively extended structures. Figure [1] (left) shows the 1.3-mm continuum map. It is seen that 
MM1 and MM2 are resolved into at least three and two smaller mm peaks res pectively (marked 



by cro sses in Fig. [I]). Only MM1 is associated with the cm emission detect ed by|R engaraj an &: Ho 



1996). Two of the three faint 7-mm sources (Sj mm ~ 1 mJy) reported by Ivan der Tak &: Menten 
(|2005l ) toward MM1 at a resolution of ~ 0.05" (marked by triangles in Fig. Q]) are counterparts 
of the 1.3-mm peaks. The faintest 7-mm source has no association in our continuum or line data. 
In addition to the clearly identified 1.3-mm peaks, the northeast-southwest large-scale continuum 
ridge appears to have more fainter sources. Another possible source is well separated from the 
ridge, at ~ 8" to the southwest of MM1. More sensitive observations are necessary to investigate 
their nature. 

We label the identified mm peaks as MMl-Main (the brightest source of MM1), MM1-NW (for 
northwest), MM1-SE (for southeast), MM2-Main, and MM2-NE. Table 2 lists the peak positions 
and peak intensities measured in the mm map of Fig. Q] left. The sums of the 1.3-mm fluxes of the 
components that we obtain from multi-component Gaussian fits to the sources comprising MM1 
and MM2 are robust, and consistent with the fluxes measured by integrating the intensity over 
the areas of interest. However, the sizes and fluxes of the individual components in the fits are 
not accurate, mainly because of insufficient angular resolution. Table 2 lists the added flux of the 
subcomponents of MM1 and MM2. The ratio of the 1.3-mm flux of MM1 to that of MM2 in our 



data is 1.2, very close to that report ed by Ivan der Tak et al.l (2000l ): 1.3. The fluxes that we report 



are 88% to 106% larger than those in Ivan der Tak et al.l (|2000l ). probably due to differences in (u, v) 
coverage and flux-scale uncertainties. 

Only the bright, compact sources are detected in the 0.9- mm continuum image (Fig. fright). 
This single-configuration data set has a more modest (u, v) coverage than the concatenated 1.3-mm 
data and is less sensitive to extended structures. 
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3.1.2. Nature of the Continuum 



To set an upper limit to the free-free contribution at 1.3-mm we extrapolate the 8.4 to 43.3- 
GHz free-free spe ctral index a = 1.03 ± . 08 (w here the flux goes as S v oc v a ), calculated from the 



fluxes reported by van der Tak Sz Menten ( 2005 ) and Rengaraian Hq (1996). This i s a reasonable 
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For MM1 with a 7-mm flux of Sj mm ~ 4.3 mJy the maximum free-free flux at 1.3-mm is 28 
mJy The 1.3-mm flux integrated over MM1 is Si. 3mm ~ 357 ± 25 mJy (Table 2). Therefore, the 
free-free contribution to the 1.3-mm flux is at most ~ 8 %. No cm continuum has been detected 
toward MM2 or in the rest of the field, thus, the 1.3-mm emission outside MM1 is most probably 
produced entirely by dust. Using the same considerations, the free-free flux of MM1 at 0.9 mm 
is less than 42 mJy. The integrated flux of MM1 at this wavelength is 5*0. g m m ~ 612 mJy, then 
the free-free emission is at most y 7 %. Because the data are taken at different epochs, the 
possibility of radio variabi lity (see iFranco-Hernandez &: Rodriguez] 12004 : Ivan der Tak et al.1 120051 : 
Galvan-Madrid et al.ll2008l . for reports on other targets) adds to the uncertainty. The 0.9-mm data 
may suffer from missing flux, making the fractional free-free contribution at this wavelength even 
smaller. In the rest of the paper, we consider the (sub)mm free-free emission to be negligible. 

To compare the 0.9-mm and 1.3-mm fluxes in a consistent way, we produced images with 
a uniform (u, v) coverage (using only baselines with lengths from 30 to 230 KX), without self- 
calibration (using only the VEX data at 1.3 mm and the Extended data at 0.9 mm), and a common 
circular synthesized beam (HPBW = 0.85"). The average spectral indices of the two main mm 
cores are («mmi) = 3.3 ± 0.3 and («mm2) = 2.5 ± 0.4. In the Rayleigh-Jeans (R-J) approximation 
(hu <C ksT), the spectral index of thermal dust emission is a = 2 + /?, where (3 is the exponent of 
the dust absorption coefficient. The fiducial interstellar-medium (ISM) value of j3 is 2, while for 
hot cores in massive star-forming regions (MSFRs) typical values are j3 ~ 1 — 2 (e.g., Cesaroni et 
al. 1999, Zhang et al. 2007). Therefore, MM1 has (3 « 1.3 typical of a hot core, but MM2 has 
f3 ~ 0.5. In Section \3. 21 we show that the kinetic temperature of MM2 is ~ 46 K, then the R-J limit 
is not a good approximation at 0.9 mm for MM2. 

Without using the assumption of being in the R-J limit, the gas mass M gas derived from 
optically-thin dust emission at 1.3-mm can be obtained from Kirchhoff's law: 



M, 



Mr, 



(26.6) x exp 



11.1 



PWK] 



[Fi .3mm /Jy][d/kpc]< 

[«1.3mm/cm 2 g _1 ] 



(1) 



where 7d ust is the dust temperature, Fi^ mrn is the 1.3-mm flux density, d is the distance to the 
object, and Ki.3mm is the dust absorption coefficient. Assuming coupling between gas and dust, the 
dust temperature in MM1 is ~ 347 K in the inner arcsecond (obtained from fits to CH3CN lines, 
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see Section [3.2.2p and > 100 K at larger scales (obtained from NH3 lines, see Section 13.2.11) . Fo r 
this range of temperature, using an opacity Ki,3 mm = 0.5 cm 2 g _1 (jOssenkopf Henning]|l994l ). 
Equation (1) gives a mass for MM1 in the range Mmmi = [9, 32] Mq. For MM2 with a temperature 
of 46 K (Section I3.2.ip . the mass is Mmm2 ~ 60 Mq. MM2 then appears to be much colder and 
more massive (in gas) than MM1. The uncertainties in opacity make the mass estimation accurate 
to only within a factor of a few. 



3.2. Molecular Line Emission 

3.2.1. The Par sec- scale Gas 

The large-scale gas within an area of ~ 1' x V (or ~ 1 pc) can be divided into quiescent gas 
and high-velocity gas. The quiescent gas is best traced by the VLA NH3 data. The high-velocity 
gas is seen in the SMA CO (2-1) maps. 



Morphology and Velocity Structure 



There is a clear morphological difference between the quiescent and the high- velocity gas. Figure 
[2al shows moment maps of the NH3 (2,2) line overlaid with the high- velocity CO gas. The NH3 
moment maps were integrated in the [31,43] km s _1 LSR velocity range. The blueshifted CO 
gas was integrated in the range [0,22] km s~\ and the redshifted CO was integrated in [62, 98] 

km s _1 . The systemic velocity of the gas closer to MM1 is V sys ~ 38.5 km s _1 (Section I3.2.2[) . 
The quiescent NH3 emission is composed of one prominent filamentary structure in the east-west 
direction that peaks toward the MM1 region (Fig. I2a|) . plus another filamentary structure that 
extends to the south of MM1 and MM2, and some fainter clumps toward the northwest of MM2. 
The high-velocity CO traces at least two molecular outflows that expand outward off the quiescent 
filaments. The lobes of the most prominent outflow are centered in MM1, and extend toward 
the northwest (redshifted gas) and southeast (blueshifted) at a position angle P. A. ~ 133°. The 
observed size of this outflow is about 0.4 pc. The redshifted lobe of a second high-velocity outflow 
extends « 0.5 pc to the north-northeast of the cores at P. A. « 19°, and appears to be originated in 
MM2. The blueshifted side of this outflow does not appear at high velocities. The P.A. that we find 
for the main outflow a grees very well with the P.A. ~ 135° of the outflow as seen at 2.2 ^m reported 
by (|Davies et al.ll2010l. see their Figure 1), w hich matches an elongated 4.5 fim structure in Spitzer 
images (Figure 12 of lCvganowski et al.ll2008P). However, the infr ared emission is three to four times 
larger. Also, using single-dish observations, Ide Wit et al.l (|2010i ) reported a CO J = 3 — 2 outflow 
whose orientation matches those of both the SMA and the near IR outflows. The low-velocity CO 
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(2-1) gas could not be properly imaged because of the lack of short (u, v) spacingqj. 

The large-scale quiescent gas has two velocity components, which appear to be two different 
structures of gas at different velocities, as can be seen in Fig. 2b. The gas associated with MM1 
and MM2 (center of the main filament), as well as the western part of the main filament and the 
north-south extensions appears to be at ~ 38.5 km s _1 , with a typical mean- velocity dispersion 
of 0.4 km s _1 . The eastern part of the main filament appears blueshifted, with a mean centroid 
velocity of ~ 35.9 km s _1 and a dispersion of the centroid velocity of 0.4 km s" 1 . The two gas 
structures overlap in space toward the MM1/MM2 region, which suggests that the star formation 
activity in these mm cores was triggered by the convergence of the filaments of molecular gas. Fig. 
[3] shows a three-dimensional rendering of the same data that better illustrates this result. At the 
position of the mm cores, the filamentary structures do not merely superpose in position-position- 
velocity space, but merge into a region that suddenly extends to higher velocities. The larger 
velocity range at the center (see also Fig. 2c) is found to be due to coherent velocity structures 
(disk and outflows in the dense gas) with the SMA data (Section I3.2.2H . 



Physical Parameters 



Now we derive the temperature structure of the parsec-scale filaments and lower limits to the 
outflow parameters. 

Most of the gas in the pc-scale filaments, including the gas associated with MM2, is cold, with 
a kinetic temperature Tk; n = [20,50] K. T^ n rises significantly only toward MM1. A temperature 
map at the resolution of the NH3 data ( Fig. B l) was obtained by fitting the (1,1) and (2,2) line 
profiles as described in Rosolowskv et al. ( 20081 ). The errors in the temperature determination are 
in general ~ 3 K, but get too large toward MM 10. We determine the temperature of MM2 to be 
Tmm2 ~ 46 K. The temperature of MM1 is constrained to Tmmi > 100 K. 



For the CO (2-1) line, the interferometric data suffer from missing flux for the more extended 
emission close to V^ ys . We set the following limits to the outflow parameters: mass M out > 27 M@, 
momentum P out > 233 M & km s _1 , and kinetic energy -E ut > 3 x 10 46 erg s -1 , where we corrected 
for the optical depth at each velocity bin using the 13 CO (2-1) line. We refer the reader to 
Qiu et al.l (|2009i ) for a description of the procedure to calculate the aforementioned quantities. 
The momentum and especially the energy estimations are less affected by missing flux since they 
depend more on the high- velocity channels. Ide Wit et al.l (|2010l ) estimated the inclination angle of 



3 Single-dish data were taken for this purpose but these were corrupted due to a bad off position and could not be 
used. 

4 The (2,2) to (1,1) ratio is not sensitive to temperatures much larger than 50 K, but we confirm the large 
temperatures in MM1 at smaller scales using the CH3CN lines, see Section 13.2.21 The errors in the fits also increase 
toward MM1 due to its wider velocity structure. 
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the inner-cavity walls of the outflow to be i ~ 50°, by radiative-transfer modeling of the mid-IR 
emission. Correcting by inclination, the lower limits to the momentum and energy of the outflow 
are P out > 362 M km s" 1 and E out > 7 x 10 46 erg s _1 . 

3.2.2. The Inner 0.1 pc 
Morphology and Velocity Structure 

The SMA data permit us to study the molecular gas at a resolution of « 1500 AU (0.4"). 
Some molecular lines trace relatively cold gas, while some other lines trace the warmer gas closer 
to the heating sources. Figure [5] shows the spectra over the entire sidebands from the pixels at the 
positions of the mm peaks MMl-Main and MM2-Main. The prominent lines are labeled, and listed 
along with their upper-level energy in Table 3. Lines with a peak intensity below 20 K are not 
listed. A complete inventory of the molecular lines in W33A will be presented in the future. It is 
immediately seen that MM1 has a "hot-core" spectrum, while MM2 is almost devoid of molecular 
emission, if not for the CO, 13 CO and C 18 J = 2 - 1, and faint SO J(K) = 6(5) - 5(4) emission. 
We interpret this difference as a signature of the evolutionary stage of the cores, MM1 being more 
evolved than MM2. 

Figure [6] shows the moment maps for three lines that exemplify what was mentioned above. 
The SO 6(5)-5(4) line (top row) extends in a ridge of « 0.1 pc long in the northeast-southwest 
direction, from MM1 to MM2. The emission is stronger toward MM1 and peaks in MMl-Main. 
Part of the emission toward MM2 is redshifted by ~ 1 — 2 km s _1 with respect to the emission 
in the MM1 side, but there is no clear velocity pattern. Lines such as SO likely have large optical 
depths and trace only the surface of the emitting region, where clear velocity gradients, especially 
of rotation, may not be expected. From the SO data we constrain any velocity difference between 
the MM1 and MM2 cores to AV < 2 km s~ x . 

For a given molecule, the isotopologue lines and the lines with upper energy levels well above 
100 K trace the more compact gas, closer to the heating sources. Figure [6] shows the examples of 
!3CS J = 5 - 4 (middle panel) and CH 3 CN J(K) = 12(3) - 11(3) (bottom). Both of them are only 
visible toward MM1, and peak in MMl-Main. These lines trace a clear velocity gradient centered on 
MMl-Main, the blueshifted emission is toward the southwest, and the redshifted emission toward 
the northeast, perpendicular to the main bipolar CO outflow. We interpret this as rotation. The 
emission tracing this velocity gradient is not isolated, there is also emission coming from MM1- 
NW and MM1-SE. Especially in the CH3CN lines, this extra emission appears to trace redshifted 
and blueshifted emission respectively. One possibility is that MM1-NW and MM2-SE are separate 
protostars from MMl-Main. However, the orientation of the lobes in the high-velocity outflow 
is the same. Therefore, we prefer the interpretation that MM1-NW and MM2-SE are not of 
protostellar nature, but emission enhancements (both in continuum and line emission) from the 
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hot base of the powerful molecular outflow driven by the dis k-like structure surroun d ing M M1- 
Main. In this scenario, the other 7-mm sources reported by Ivan der Tak &: Mentenl ((2005) (or 
at least the counterpart of MM1-NW) can be interpreted as shocked free-free enhancements in a 
protostellar jet, similar to those observed in the high-m ass star formation region IRAS 16547 — 4247 
((Rodriguez et all 1200a : Franco-Hernandez et al.ll2009l ). 



Figure [7] shows the position-velocity (PV) diagrams of the CH3CN K = 3 line shown in Fig. [6j 
centered at the position of MMl-Main perpendicular to the rotation axis {top frame) and along it 
{bottom frame). The rotation pattern is similar to those observed in objects that have been claimed 
to be Keplerian disks, i.e., structures where the mass of the central object is large compared to 
the mass of the gas, rotat ing with a velocity V ro t oc r" -°- 5 (|zhang et al.lll998l : Icesaroni et alJbood : 
Jimenez-Serra et al.ll2009i ). The large velocity dispersion closest to MMl-Main (Figs. [6] and [7]) 
ought to be caused by unresolved motions in the inner disk, since velocity dispersions well above 1 
km s _1 cannot be due to the gas temperature. 



Recently, iDavies et al.l (|2010i ) reported a possible disk-jet system centered in W33A MMl- 
Main. The jet, observed in the Br7 line, extends up to ±300 km s _1 in velocity at scales ~ 
1 AU, with a similar orientation and dir ection to the molec ular outflow reported in this paper. 
However, the velocity structure of what IDavies et al.l (|2010l ) interpret as a rotating disk has a 
similar orientation but opposite sense of rotation as the disk that we report. They used CO 
absorption lines with upper energy levels E u ~ 30 K, while we use emission lines like those of 
CH3CN, with E u > 100 K (Table 3). If an extended screen of cold gas with a negligible velocity 
gradient is between the observer and the inner warm gas with a velocity gradient, it is possible 
that the absorption lines are partially filled with emission, mimicking a velocity gradient with the 
opposite sense than that seen in the emission lines. 



Physical Parameters 

Now we derive the temperature and column density of the hot-core emission, and constrain 
the stellar mass, gas mass, and CH3CN abundance in MM1. 

The kinetic temperature of the innermost gas can be obtained from the K lines of CH3CN 
J = 12 — 11. To avoid the simplification of considering optically-thin emission assumed in a 
population-diagram analysis, we fit all the K lines assuming LTE, while simultaneously solving for 
the temperature lk in , column density of CH3CN molecules Ach 3 cn> an d line width at hal f-power 



FWHM. The procedure to obtain the level populations can be found in lArava et al.l (|2005l ) 



Figure [8] shows the results of our fits to the CH3CN spectra. The systemic velocity V sys fa 
38.5 km s _1 was found to be optimal. The data outside the lines of interest have been suppressed 
for clarity, and the fit was done in the frequency windows where only the lines of interest are 
present. The gas is warmer (by 18%) and denser (by 415% in column) toward the peak MMl-Main 
{bottom frame) than in the average of the sources composing MM1 {top frame). This makes the 
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case for the internal heating and a centrally-peaked density gradient in MMl-Main, as well as for 
its protostellar nature. Some lines are not completely well fit under the assumption of a single 
value for the parameters. On the one hand, the bright K = 7 and 8 lines indicate the presence of 
some column of very warm gas; on the other hand, the K = 3 is not brighter than the K < 3 lines 
within the uncertainties, indicating that some column of gas is below 200 K. Two-component fits 
do not give better results. The reported values should be interpreted as an average along the line 
of sight. Detailed radiative transfer in the context of a physical model is currently under way. 

Under the interpretation of edge-on rotation as the cause of the velocity gradient observed in 
Figs. [6] and [71 the dynamical mass in MMl-Main (stars plus gas) is about 9 M Q , where a mean- 
velocity offset of 2.0 km s -1 with respect to V^ ys ~ 38.5 km s _1 was taken at a radius of 0.5". 
Assuming that the disk is perpendicular to the outflow, and correcting for the outflow inclination 



angle i = 50° estimated bv lde Wit et al.1 (120101 ) . the enclosed dynamical mass is M<j yn ~ 15 Mq. 



Given that W33A is fragmented into multiple sources, a strict upper limit to the stellar mass 
in MMl-Main is < 20 Mq, the mass necessary to account for the total luminosity (~ 10 5 Lq) 
of W33A. Also, the total gas mass in the MM1 sources inferred from dust emission is ~ 10 Mq, 
therefore the gas mass in the rotating structure around MMl-Main should be a fraction of it. 
For a gas mass in the rotating structure M gas ~ 5 Mq, the mass of the protostar(s) in MMl-Main 
amounts to M* ~ 10 Mq. MMl-Main appears to be B-type protostar still accreting from a rotating 
disk-like structure. 

From the ~ [9,32]M of gas mass in MM1 we derive an average column density of molecular 
gas of about [1,4] x 10 23 cm -2 . For an average CH3CN column of ~ 8 x 10 15 cm -2 , a CH3CN 
abundance with respect to H2 of Xf CH.sCN) ~ [2,8] x 10~ 8 is derived, sim i lar estimates in other 



regio ns of high-mass star formation (jWilner et al.lll994l ; iRemiian et al.ll2004l ; iGalvan-Madrid et al 
2009|). 



Discussion 



4.1. Star Formation from Converging Filaments 

In W33A, two localized regions of star formation (MM1 and MM2) separated by 0.1 pc are 
surrounded by a common filamentary structure of ~ 1 pc in length (Section 13.2.1)) . The two 
velocity components of this filamentary structure are separated by ~ 2.6 km s -1 in line-of-sight 
velocity and intersect in projection right at the position of the star formation activity (Fig. 2). 
The velocity components are not a mere superposition in position-position-velocity space, but they 
merge into a structure with larger motions (Fig. [3j Section I3.2.ip . resolved into a disk/outflow 
system by subarcsecond resolution observations (Figs. [6] and [7J Section 13. 2. 2ft . This suggests 
that star formation in W33A was triggered by the con vergence of molecula r filaments. Such a 
scenario has been suggested for the region W3 IRS 5 by iRodon et al.1 (|2008l ) . More recently, the 
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merging of filaments has also been claimed bylJimenez-Serra et al.l (|2010l ) in the infrared dark cloud 
G35.39-0.33 and bv lCarrasco-Gonzalez et al.l faoid ) in the MSFR W75. 



This mode of st ar formation is predicted by numerical simulations of star formation triggered 



by co nverging flows ( Ba llesteros-Paredes et al.lll999l : iHeitsch et al.ll2008l : IVazquez-Semadeni et al 



20071 ). In those simulations, the formation of molecular clouds itself is a product of the convergence 
of streams of neutral gas. Later in the evolution of the molecular clouds, filaments of molecular gas 
can converge (merge) with each other, leading to the formation of cores and stars. 

We pres ent here a simple compar is on wi th a region found in the numerical simulation recently 
reported by IVazquez-Semadeni et al.1 (120091 ) . This simulation represents the formation of a gi- 
ant molecular complex from the convergence of two streams of warm neutral gas, at the scale of 
tens of parsecs. Specifically, th e simulation was per formed using the smoothed particle hydrody- 
namics (SPH) code GADGET (ISpringel et al.l 120011). including si nk-particle and radiative cooling 



prescriptions (jJappsen et al.l 120051 ; IVazquez-Semadeni et al.l 120071 ) . The convergence of the warm 



diffuse-gas streams triggers a thermal instability in the gas, which causes it to undergo a transition 
to the cold atomic phase, forming a cloud. The latter soon becomes gravitationally unstable, begins 
contracting, and undergoes hierarchical fragmentation. During the contraction, the density of iso- 
lated clumps increases and they can reach physical conditions corresponding to those of molecular 
clumps. Finally, the global collapse reaches the center of mass of the cloud, at which point a region 
with physical conditions corresponding to those of MSFRs forms. The simulation box has 256 pc 
per side, and the converging flows have a length of 112 pc, and a radius of 32 pc. However, after 
the gravitational contraction, the clumps are only a few parsecs across. Since SPH is essentially 
a Lagrangian method, it allows suffic ient resolution in these dense regions. We refer the reader to 
Vazquez-Semadeni et al.l (120071. 12009!) for det a ils of the simulation. Here we focus on the region 



Vazquez-Semadeni et al.l (|2009i ). albeit roughly 1.5 Myr later than the time 



called "Cloud 1" in 
examined in that paper. 

Figure [9] shows two snapshots of column density separated by 0.133 Myr (the time interval 
between successive data dumps of the simulation). The column density is computed by integrating 
the density along the x-direction over the 10 pc path 123 < x < 133 pc, which is centered at the 
midplane of the simulation, where the (sheet-like) cloud is located. In this region and epoch of the 
simulation, the two leftmost filaments in the top panel of Fig. [9] converge to form the boomerang- 
shaped filament seen in the bottom panel. Note that the simulation was not designed to simulate 
the observed filaments. The comparison is only intended to show that some of observed properties 
of W33A can arise naturally in the context of a simulation of the formation of a large molecular 
complex that contains a filament system. 

The peaks of two filaments of gas are initially separated by 0.4 pc {top frame of Fig. [9]) and 
then merge at a projected velocity of 3 ± 1.5 km s , measured directly from the displacement 
observed between the two panels. The column density of the filaments is in the range [10 3 , 10 4 ] code 
units (Fig. [9]), or Nr 2 = [0-5 x 10 24 ,0.5 x 10 25 ] cnr 2 . For an NH3 abundance with respect to H2 in 
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the range X(NH 3 ) = [1(T 8 , 10~ 7 ] (e.g., iGalvan-Madrid et al.ll2009h . the column density of the cold 
filaments in our observations is Nr 2 = [1-1 x 10 24 ,1.1 x 10 25 ] cm -2 . Figure PTOl shows the volume 
density (color scale) and y— z plane velocity (arrows) in a slice trough the filaments. It is seen 
that the filaments reach densities typical of MSFRs (peak n ~ 10 5 cm -3 ) and that their velocity 
field presents fast jumps of a few km s _1 in the interaction zones, comparable to our observations. 
We conclude that some of the properties of the observed filaments such as sizes, column densities, 
and velocities agree within a factor of 2 with those from the simulation. This rough comparison 
illustrates that our interpretation of convergence between the observed filaments is feasible. 



4.2. Cores at Different Evolutionary Stages 



The star-forming cores in W33A appear to be at markedly different evolutionary stages (Section 
I3.2.2p . The first piece of evidence for this is the clear difference in the richness of the molecular- line 
emission from MM 1 to MM2 (see Fig. [5]). MM1 has molecular emission typical of a "hot-core" 
(jKurtz et al.l 12000 ) . with a prominent CH 3 CN J = 12 — 11 "if -forest" that can be detected up to 
the K = 8 line, with upper-level energy E u = 525.5 K (see Table 3). The average gas temperature 
of MM1 is ~ 347 K. In contrast, MM2 is almost devoid of "hot-core" emission, and is only detected 
in a few molecules. MM2 is much colder than MM1, with a temperature Tmm2 ~ 46 K. The 
second piece of evidence is the mass content of the cores. MM1 has only ~ [9, 32] Mq of gas, while 
MM2 has ~ 60 Mq. This could naively be interpreted as MM2 having a much larger gas reservoir 
than MM1, but it should be kept in mind that both cores are part of a common parsec- scale 
structure. Clu mp infall at pc scales has be en reported from single-dish (|Wu Evansl 120031 ) and 
interferometric (jGalvan-Madrid et al.ll2009l ) observations. Also, numerical simulations of parsec- 
sized clumps show that the star-forming cores that give b irth to massive sta r s are continuously 



fed from gas in the environment at the clump sc ales (e.g., iBonnell et al J 12003 ; ISmith et al 



2009 



Vazquez-Semadeni et al. 



2009 



Peters et alJl201(l ). 



There are two possibilities that we briefly discuss here: (1) the prestellar cores that were the 
precursors to MM1 and MM2 appeared at the same time and then MM1 evolved faster to produce 
a ~ 10 Mq star (Section 13.2.2)) . while MM2 only produced at most an intermediate mass star 
(MM2 is not prestellar, since it has some internal heating and appears to power an outflow, see 
Section I3.2.ip . or (2) the prestellar core precursor to MM2 formed later and has yet to form at 
least one massive star. Our observations cannot tell these options apart. A measurement of the 
accretion rate in both cores would be helpful. Sources at different evolution ary stages within a 
single star- for ming cluster have also been reporte d recently for AFGL 5142 (jZhang et all 120071 ). 
G28.34+0.06 rtzhang et alJboOfll ). and AFGL 961 (|williams et alJboOfll ). 
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4.3. A Rotating Disk/Outflow System in MMl-Main 



In the past decade, the question of whether massive stars form by disk-outflow mediated 
accretion similar to low-mass stars has been the subject of intensive research. The answer is 
positive: they definitely do. Some of the massive protost ars that have been shown to harbo r 
disk/outflow system s are G192. 16-3 . 82 dsheoherd et all hoOlh. Cep heus A HW2 rtPatel et al lbpOfih . 



IRAS 20126+4104 (|Cesaroni et alJboOfil ). and IRAS 16547-4247 (|Franco-Hernandez et al 



2009) 



All these relatively clean disk examples, however, do not have stars more massive than 15 to 
20 Mq. More massive (O-type) stars have also been shown to form via disk- mediated accretion. 
However, the innermost part of the accretion flow is often (at least partially) ionized and is observed 
as a hypercompact (HC) H II region. Also, the gas in these more massive regions is warmed 
up to farther distances, and very massive rotating st ructures of size up to 0.1 pc are typically 
observed. Examples a re G10. 6-0.4 dKeto Woodlbood'l . G2 4.78+0 08 feeltran et alJ boodl. W51e2 
(|Klaassen et all l2009h . and G20.08-0.14 N dGalvan-Madrid et all 120091 ) . all of whic h have stellar 
masses above 20 Mq. A possible exception to this scenario is W51 North, where IZapata et al 



(120081 . l2009i ) claims to have found a protostar with M* > 60 Mq and without a "bright" (with 
flux above tens of mJy at wavelengths ~ 1 cm) H II region. This apparent discrepancy is solved if 
multiple, lower-mass stars account for the mass in W51 North or if the H II region in this source is 
gravitationally trapped as currently observed. Indeed, detailed simulations of the evolution of HC 
H II regions within accretion fl ows show that thei r radio-continuum emission flickers significantly 
in timescales from 10 to 10 4 yr ( Peters et al.ll2010l ). 



In this paper we report on the existence of a rotating disk centered on MMl-Main in W33A 
(Section 13.2.2)) . To our sensitivity, the radius of the disk is R < 4000 AU. The warm-gas emission 
does not come only from the disk, but also from an structure elongated perpendicular to it, coin- 
cident with the other mm peaks (MM1-NW and MM1-SE) along the direction of the outflow (see 
Figs. [6]and[7|. We propose that the secondary mm peaks in MM1 are not of protostellar nature, 
but regions where the emission is enhanced due to the interactio n of the outflow with the disk and 
its inner envelope. Indeed, the two brightest 7-mm detections of Ivan der Tak fe Mentenl (|2005l ) are 
counterparts of MMl-Main (Ql, S v « 1.7 mJy at 7 mm) and MM1-SE (Q2, S v ps 0.6 mJy at 7 



mm . 



The free-free emission from the 7 mm source Ql should be a combination of photoioniza- 
tion by the central protostar and shock-induced ionization of material due to the jet observed by 
Davies et al .1 ()2010l ). likely dominated by the latter. The origin of the free-free emission from Q2 
could also be shocks, although deeper 7 mm observations are needed to confirm this source. Q3 
does not have a mm/submm counterpart and may not be a real s ource. The radio con t inuum 



sources of W33A could then be analogs of those in IRAS 16547-4247 ((Rodriguez et al 



2005 



2008) 



which have fluxes a factor of a few larger at 0.75 times the distance to W33A. Indeed, we find that 
the radio-continuum emission from Ql (MMl-Main) agrees with the correlation found for low-mass 
jets between the radio-continuum lumin osity of the jet and the momentum rate of the associated 
molecular outflow: P = 10 _2 ' 5 (S , iy (i 2 ) 1 ' 1 ( Anglada et al.lll998l ). where the radio luminosity has units 
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For W33A with a radio flux 



of mJy kpc 2 and t he momentum rate has u nits of M@ yr _1 km s~ 
<S3.6cm = 0.79 mJy (jRengaraian &: Holll996l ). and a distance d = 3.8 kpc (jJaffe et al.lll982l ). the ex- 
pected momentum rate is P = 0.046 M Q yr -1 km s _1 , while the observed momentum rate (lower 
limit) is P = 0.040 M yr _1 km s _1 (obtained dividing the momentum of the molecular outflow by 
its length, including the inclination correction). This correlation was found to ho ld for three well 



studie d massive protostars (IRAS 16547-4247, HH 80-81, and Cep A HW2) bv iRodrfguez et al 



(|2008l ). In this paper we report that it also holds for W33A, which constitutes further evidence for 
a common accretion mechanism between low- and high-mass protostars, at least to the stage prior 
to the development of a brighter H II region. 

Comparing the ratio of the radio luminosity to the IR luminosity L(8 GHz)/L(IR) with the 
recombination-line line width may be a useful criterion to distinguish between a source ionized 
by shocks (jet or stellar wind) or by photo ionization (what usually is called an H II region). In 
Fig. 6 of the review bv iHoare et al.l (|2007l ). it is seen that UC H II regions have the largest L(8 
GHz)/L(IR) and the smallest line width, jet sources have the smallest L(8 GHz)/L(IR) and the 
largest line width, and HC H II regions fall in between the previous two. For W33A, log(L(8 
GHz)/L(IR)) ~ 7.1, and the FWHM of the IR recombination lines is of several hundreds km s _1 , 
again consistent with a jet source. 

The star(s) at the center of MMl-Main (with ~ 10 M ) appears to dominate the dynamics 
of the disk, but we cannot rule out the existence of additional, less massive objects within it. Indeed, 
models of massive protostel lar disks predict their fragmentation and the fo r mation of a few lo wer 



mass companions within it (jKratter &; Matznerl 12009: 



Krumholz et al. 



2009 



Peters et alJl201(J ). 



5. Summary and Conclusions 

We present for the first time resolved observations in both mm continuum and molecular-line 
emission for the massive star formation region W33A, characterized by a very high luminosity 
[L ~ 10 5 Lq) and very low radio-continuum emission (~ 1 mJy). Both of the previously known 
mm cores (MM1 and MM2) are resolved into multiple peaks, and appear to be at very different 
evolutionary stages, as indicated by their molecular spectra, masses, temperatures, and continuum 
spectral indices. The brightest core (MMl-Main at the center of MM1) is centered on a very faint 
free-free source and the gas dynamics up to a few thousand AU of it indicates the presence of 
a circumstellar disk rotating around a stellar mass of M* ~ 10 Mq. MMl-Main also drives a 
powerful, high-velocity molecular outflow perpendicular to the disk. MM2, the coldest and most 
massive core, is not detected in hot-core lines but appears to drive a more modest outflow. Both 
MM1 and MM2 are located at the intersection of parsec-scale filamentary structures with line- 
of-sight velocity offset by ~ 2.6 km s -1 . Analysis of the position-position-velocity structure of 
these filaments and a comparison with recent numerical simulations suggests that star formation 
in W33A was triggered by the convergence of filaments of cold molecular gas. 
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Table 1. Observational Parameters 



Epoch 


Array 


Phase Center a 


Bandpass 


Flux 


Phase 






a(J2000) 


5(J2000) 


Calibrator 


Calibrator 


Calibrator 


2004 06 14+15 


VLA-D 


18 14 39.500 


-17 51 59.800 


3C273 


3C286 


1851 + 005 


2006 05 22 


SMA-Extended 


18 14 39.509 


-17 51 59.999 


3C273 


Callisto 


1733 - 130 


2007 07 17 


SMA-Compact 


18 14 39.495 


-17 51 59.800 


3C454.3 


Callisto 


1733 - 130 


2008 08 02 


SMA-VEX 


18 14 39.495 


-17 51 59.800 


3C454.3 


Callisto 


1733 - 130 



a Units of RA. are hours, minutes, and seconds. Units of decl. are degrees, arcminutes, and arcseconds. 



Table 2. Millimeter Continuum Sources 



Core a 


Component 15 


a(J2000) c 


<5(J2000) C 


-fpoak(l-3 mm) d 


5(1.3 mm)° 






(hrs,min,sec) 


(deg,arcmin,arcsec) 


( mjy beam -1 ) 


(mjy) 




MM1-NW 


18 14 39.47 


-17 51 59.7 


31 


357 


MM1 


MMl-Main 


18 14 39.51 


-17 52 00.0 


65 


357 




MM1-SE 


18 14 39.55 


-17 52 00.4 


25 


357 


MM2 


MM2-Main 


18 14 39.24 


-17 52 01.9 


43 


289 




MM2-NE 


18 14 39.31 


-17 52 00.6 


22 


289 



a Main core as labeled in Fig. [T] left panel. 

b Clearly distinct subcomponents of the main cores as marked in Fig. [T] 
c Position of peak. 

d Peak intensity ±1.5 mjy beam" 1 . HPBW = 0.63" X 0.43". 

c Added flux of the subcomponents of each core. The uncertainties in the fluxes of MM1 and 
MM2 are ±20 mjy and ±25 mjy respectively. 
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Table 3. Bright Molecular Lines' 



Species 


Transition 




E-u 






(GHz) 


(K) 


c 18 o 


2-1 


219.5603 


15.8 


HNCO 


10(2,9)-9(2,8) 


219.7338 


228.4 


- 


10(2,8)-9(2,7) 


219.7371 


228.2 


- 


10(0,10)-9(0,9) 


219.7982 


58.0 


H2 13 CO 


3(1,2)-2(1,1) 


219.9084 


32.9 


SO 


6(5)-5(4) 


219.9494 


34.9 


CH3OH 


8(0,8)-7(l,6) E 


220.0784 


96.6 


13 co 


2-1 


220.3986 


15.86 


CH3CN 


12(8)-11(8) 


220.4758 


525.5 




12(7)-11(7) 


220.5393 


418.6 


HNCO 


10(1,9)-9(1,8) 


220.5847 


101.5 


CH 3 CN 


12(6)-11(6) 


220.5944 


325.8 




12(5)-11(5) 


230.6410 


247.3 




12(4)-11(4) 


220.6792 


183.1 




12(3)-11(3) 


220.7090 


133.1 




12(2)-11(2) 


220.7302 


97.4 




12(1)-11(1) 


220.7430 


76.0 




12(0)-11(0) 


220.7472 


68.8 


CH3OH 


15(4,11)-16(3,13) E 


229.5890 


374.4 




8(-l,8)-7(0,7) E 


229.7588 


89.1 




3(-2,2)-4(-l,4) E 


230.0270 


39.8 


12 CO 


2-1 


230.5380 


16.5 


ocs 


19-18 


231.0609 


110.8 


13 cs 


5-4 


231.2207 


26.6 


CH30H 


10(2,9)-9(3,6) A— 


231.2811 


165.3 



a Molecular lines with peak T B > 20 K. The first 
column refers to the molecule tag, the second col- 
umn to the transition, the third to its rest frequency 
as found in Splatalogue (http:/ /www. splatalogue.net/ ), 
and the fourth to the upper-level energy. Data used 
from Splatalogue a re compiled from the C DMS catalog 
jMiiller et alj|2005ft and the NIST catalog llLovaslbOoA 
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Fig. 1. — (Sub) millimeter continuum emission in W33A. The left panel shows the 231 GHz (1.3 
mm) continuum from the SMA Compact +VEX data (HPBW = 0.63" x 0.43", P.A. = 30.7°). 
Contours are at —5,5,7,10,15,20,30, and 40 times the noise of 1.5 mJy beam -1 . The right 
panel shows the 336 GHz (0.9 mm) continuum from the Extended-configuration data (HPBW 
= 0.88" x 0.83", P.A. = 275.1°), with contours at -5, 5, 7, 10, 15, 20, 30, and 39 times the rms noise 
of 6 mJy beam -1 . The cores MM1 and MM2 are labeled, and the sources into which they fragment 
are marked by crosses. Triangles mark the positions of the faint 7 mm sources reported by van der 
Tak & Menten (2005). 1 arcsec corresponds to 3800 AU (0.018 pc). 
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Fig. 2a. — Parsec-scale gas structure toward W33A. The gray scale shows the integrated intensity 
(moment 0, in mJy beam" _1 km s _1 ) of the quiescent NH 3 (2,2) gas (HPBW = 6.0" x 2.6", P.A. 
= 2°). The blue solid contours show the high-velocity gas as detected in CO (2-1) (HPBW= 3.0" x 
2.0", P.A. = 56°) integrated in the range [0,22] km s _1 . The red dashed contours show the high- 
velocity CO gas in the range [62,98] km s -1 . Contour levels are -5, 5, 7, 10, 15, 20, 25, 30, 35, 40, 50, 
and 60 x 0.7 Jy beam -1 km s _1 . The mm continuum sources identified in this paper are marked 
with stars. The directions of the identified outflows are marked with arrows. 10 arcsec corresponds 
to 38,000 AU (0.184 pc). 
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Fig. 2b (continued). — Parsec-scale gas structure toward W33A. The gray scale shows the intensity- 
weighted mean velocity (moment 1, in km s -1 with respect to the LSR) of the quiescent NH3 (2,2) 
gas (HPBW = 6.0" x 2.6", P. A. = 2°). The blue solid contours show the high- velocity gas as 
detected in CO (2-1) (HPBW = 3.0" x 2.0", P.A. = 56°) integrated in the range [0,22] km s" 1 . 
The red dashed contours show the high-velocity CO gas in the range [62,98] kms -1 . Contour 
levels are —5, 5, 7, 10, 15, 20, 25, 30, 35, 40, 50, and 60 x 0.7 Jy beam -1 km s" 1 . The mm continuum 
sources identified in this paper are marked with stars. 10 arcsec corresponds to 38,000 AU (0.184 
pc). 



-24- 




Fig. 2c (continued). — Parsec-scale gas structure toward W33A. The gray scale shows the velocity 
dispersion with respect to the mean velocity (moment 2, FWHM/2.35 in km s _1 ) of the quiescent 
NH 3 (2,2) gas (HPBW = 6.0" x 2.6", P.A. = 2°). The blue solid contours show the high-velocity 
gas as detected in CO (2-1) (HPBW = 3.0" x 2.0", P.A. = 56°) integrated in the range [0,22] 
km s _1 . The red dashed contours show the high- velocity CO gas in the range [62,98] km s _1 . 
Contour levels are -5,5,7,10,15,20,25,30,35,40,50, and 60 x 0.7 Jy beam" 1 km s _1 . The mm 
continuum sources identified in this paper are marked with stars. 10 arcsec corresponds to 38,000 
AU (0.184 pc). 
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Fig. 3. — Three-dimensional (position-position-velocity) rendering of the NH3 (2,2) data. Every 
voxel in the data cube with intensity > 20 mJy beam -1 (5a) has been included. The vertical 
axis is color-coded according to Vlsr- It is seen that the two filamentary structures at different 
velocities do not merely superpose at the position center, but merge in position-position-velocity 
space, suggesting interaction of the filaments. 
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Fig. 4. — Map of the kinetic temperature Tk; n at large scales obtained from the NH3 (1,1) and 
(2,2) data. It is seen that the pc-scale filaments are cold, with Tk; n = [20, 50] K. Only toward MM1 
Tki n rises significantly, but the errors toward this region increase up to ~ 40 K. Symbols are as in 
Figd 
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Fig. 5. — Wide-band, continuum-free spectra in the image domain from the SMA Compact +VEX 
data at 1.3 mm. The top row shows the spectra for MM1 (the peak MMl-Main). The bottom row 
shows the spectra for MM2 (the peak MM2-Main). There is a striking difference in the richness of 
the spectra between the two cores. MM2 is almost devoid of molecular-line emission, in spite of it 
having a larger gas reservoir than MM1 (Section I3.1.2p . 
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Fig. 6. — Hot-core molecules toward the center of W33A. The top row shows the SO 6(5)— 5(4) 
line. The middle row shows the 13 CS 5-4 line. The bottom row shows the CH3CN 12(3)— 11(3) line. 
Contours show the moment maps at 5, 15, 30, 50, 100, 150, and 200 x 0.05 Jy beam -1 km s _1 . The 
color scale shows the moment 1 maps (left column) and moment 2 maps (right column). Symbols 
are as in Fig. 1. While the SO traces an extended envelope reaching MM2, the other molecules 
trace a clear velocity gradient indicative of rotation centered in MMl-Main. The velocity dispersion 
also peaks in MMl-Main. 
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Fig. 7. — Position- Velocity (PV) diagrams for the CH3CN J = 12 — 11 K = 3 data. The center 
is the position of MMl-Main (Table 2). Top: PV diagram at P. A. = 39°. A clear velocity gradient 
is seen from the southwest (blueshifted, negative position) to the northeast (redshifted, positive 
position). At the angle of this cut the velocity gradient has maximum symmetry. Bottom: PV 
diagram at P. A. = 39° + 90° = 129°. Negative positions are to the northwest, and positive positions 
to the southeast. There is emission in the range [36, 41] km s _1 at all positions. Closer to the 
position center, the velocity dispersion increases rapidly. 
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Fig. 8. — Simultaneous fits to the CH3CN J = 12 — 11 K lines. The dashed gray line is the data. 
The solid black line is the fit. The data outside the lines of interest have been set to zero to avoid 
contamination by other molecular lines. Top: Average spectra over the entire MM1 region. Bottom: 
Spectra toward the few central pixels at the peak position (MMl-Main). The gas is warmer, denser, 
and has a larger line width toward this position. 
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Fig. 9. — Snapshots of molecular filaments merging with each other extracted from the simulation 
presented by Vazquez-Semadeni et al. (2009). The time interval between the first (top) and second 
(bottom) frames is 0.133 Myr. The units of the z- and y-axes are pc. The color scale shows the 
column density in code units, equivalent to 4.9 x 10 20 H2 particles per square cm. 
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Fig. 10. — Snapshots of molecular filaments merging with each other extracted from the simulation 
presented by Vazquez-Semadeni et al. (2009). The time interval between the first (top) and second 
(bottom) frames is 0.133 Myr. The units of the z- and y-axes are pc. The color scale shows the 
volume density (cm -3 ) in a slice passing through the merging filaments. The arrows indicate the 
velocity of the gas in the y-z plane. In the second frame (bottom), the two merged filaments have 
velocities roughly opposing each other by a few km s~ . 



